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We give indications that outer future trapping horizons play a role in the particular semi- 
classical instability of an evolving black hole that produces the Hawking's radiation. These 
are obtained with the use of the Hamilton- Jacobi tunneling method. It automatically selects 
one special expression for the surface gravity of a changing horizon, the one defined a decade 
ago by Hay ward using Kodama's theory of spherically symmetric gravitational fields. The 
method also applies to point masses embedded in an expanding universe and to general, 
spherically symmetric black holes. The local surface gravity solves a puzzle concerning the 
charged stringy black holes, namely that it vanishes in the extremal limit, whereas the Killing 
global gravity does not. 
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Devoted to Prof. I. Brevik on the occasion of his 70th birthday. 



1 Introduction 



It has long been felt that the usual semi-classical treatment of stationary black holes (abbrev. 
BHs) should be extended to cover at least slowly changing, or evolving black holes. By this 
expression we mean black holes that can still be described in terms of few multiplole moments 
^ such as mass, angular momentum and the charges associated to local gauge symmetries, except 

that the parameters and the causal structure change with time either because matter and gravi- 
tational radiation fall in, or because there operate a Hawking's process of quantum evaporation 
or finally because the hole is actually immersed in a slowly expanding universe. A technical 
definition of a "slowly varying BH" can be given in some cases, an example being the Booth- 
Fairhurst slowly evolving horizon, but in general it depends on the actual physical processes 
involved. For example, in the case of Hawking's evaporation, conditions for slowness in the 
presence of a near-horizon viscous fluid have been given by Brevik [1] in an interesting attempt 
to generalize 't Hooft's model of the self-screening Hawking atmosphere (quantum corrections to 
this model can be found in [2]). In general it is understood that the black hole temperature is to 
be be much smaller than the Planck mass, or equivalently the mass M » Mp = G~ 1 / 2 ~ 10 19 
Gev, while in order to study the effects of the expansion the Hubble rate H -1 should dominate 
over the black hole emission/absorption rate. 

One surprising aspect of the semi-classical results obtained so far, is that the radiation caused 
by the changing metric of the collapsing star approaches a steady outgoing flux for large times, 
implying a drastic violation of energy conservation if one neglects the back reaction of the 
quantum radiation on the structure of spacetime. But the back reaction problem has not been 
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solved yet in a satisfactory way. As pointed out by Fredenhagen and Haag long ago [3], if the 
back reaction is taken into account by letting the mass of the black hole to change with time, 
then the radiation will possibly originate from the surface of the black hole at all times after its 
formation. 

This poses the question: what is and where is the surface of a dynamical black hole? This issue 
always baffled scientists from the very beginning and produced several reactions during the 
nineties, which eventually culminated with the notion of outer trapping horizons by Hayward[4] 
and the isolated [5] and dynamical horizons of Ashtekar and Krishnan [6, 7] (a fine review is 
in [8]). Thus one is concerned to show, in the first place, what kind of a surface a dynamical 
horizon can be and also which definition can capture a useful local notion of such a surface, and 
then what sort of instability, if any, really occurs near the horizon of the changing black hole. 
This question was non trivial since a changing horizon is typically embedded in a dynamical 
space-time and it is not even expected to be a null hypersurface, although it is still one of infinite 
red shift. 

We analyze this question for a class of dynamical black hole solutions that was inspired by prob- 
lems not directly related to black hole physics, although these were subsequently reconsidered 
in the light of the black hole back reaction problem in the early Eighties. The metrics we shall 
consider are the Vaidya radiating metric [9], as revisited by J. Bardeen [10] and J. York [11], 
together with what really is a fake dynamical black holes, the McVittie solution representing, 
in author's mind, a point mass in cosmology [12]. We shall indicate how the results can be ex- 
tended to all dynamical, spherically symmetric solutions admitting a possibly dynamical future 
outer trapping horizon. 

2 Horizons 

After the time lasting textbook definition of the event horizon (abbr. EH) to be found in the 
Hawking & Ellis renowned book [13], several quasi-local notions of dynamical horizons appeared 
in the literature (a nice review is in [14]) , perhaps starting with the perfect horizons 1 of Hajicek 
[15] and the apparent horizons (AHs, boundaries of trapped 3-dimensional space-like regions) 
of Hawking-Ellis themselves. But the former only applied to equilibrium BHs and the existence 
of the latter is tie to a partial Cauchy surface so it represents only a "localization in time". 
Moreover it has proven not possible to formulate thermodynamical laws for AHs akin to those 
holding for event horizons. 

The first succesfull attempt to go beyond the limitations imposed either by the istantaneous 
character of the apparent horizons or by the global, teleological nature of the event horizons is 
due to S. Hayward. His concept of a future outer trapping horizon (FOTH) then evolved either 
into some less constrained definition, like the Ashtekar-Krishnan dynamical horizons (DH), or 
some specialization like the Booth- Fairhurst slowly evolving FOTH [16]; so an updated (but 
perhaps partial) list of locally or quasi-locally defined horizons would contain: 

(a) Trapping horizons (Hayward [4]) 

(b) Dynamical horizons (Ashtekar & Krishnan [7, 6]) 

(c) Non expanding and perfect horizons (Hajicek [15]) 

(d) Isolated and weakly isolated horizons (Ashtekar et al. [5]) 

(e) Slowly evolving horizons (Booth & Fairhurst [16]) 

1 These are null hyp ersurf aces whose rays have zero expansion and intersect space-like hypersurfaces in compact 
sets. All stationary horizons are perfect, but the converse is not true. 
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In contrast to the old fashioned apparent horizons, these newly defined horizons do not require a 
space-like hypersurface, no notion of interior and exterior and no conditions referring to infinity 
(all are non local conditions). Moreover they are not teleological and, given a solution of Einstein 
equations, one can find whether they exist by purely local computations. Finally, unlike EHs 
they are related to regions endowed by strong gravitational fields and absent in weak field regions. 
All quasi-local horizons rely on the local concept of trapped (marginally trapped) surface: this 
is a space- like closed 2-manifold S such that Om^tn) > 0; where £, n are the future-directed null 
normals to S, normalized to £ ■ n = —1, and 9^, are the respective expansion scalars. We 
write the induced metric on each S in the form 

Qab = 9ab + taUb + £bU a (2.1) 

and put q ab = g ab + £ a n b + £ b n a , not an inverse. Then q% is the projection tensor to T*(S), 
the tangent space to S. To cover BHs rather than white holes it is further assumed that both 
expansions are negative (non positive). 

The most important quantities associated with the null vector fields £ and n are the projected 
tensor fields & a b = qnQm^ah and = q^c^nS ' a n b and their decomposition into symmetric, 
anti-symmetric and trace part. Their twists are zero since they are normal to S. Finally, the 
expansions are given by 

W = q ab VJ b , 9 (n) = q ab V a n b (2.2) 

Let us describe the listed horizons in turn, adding comments where it seems appropriate. A 
black triangle down T will close the definitions. 

Future Outer Trapping Horizon: A future outer trapping horizon (FOTH) is a smooth three- 
dimensional sub-manifold H of space-time which is foliated by closed space-like two-manifolds 
St, t £ K, with future-directed null normals £ and n such that (i) the expansion 9^ of the null 
normal £ vanishes, (ii) the expansion 6^ of n is negative and (iii) C n 9^ < 0. T 
Condition (i) requires strong fields since certainly 9tn > in weak fields. Condition (ii) is related 
to the idea that H is of the future type (i. e. a BH rather than a WH); (iii) says that H is of the 
outer type, since a motion of St along n a makes it trapped. It also distinguishes BH horizons 
from cosmological ones. 

One can always found a scalar field C on H so that 

V a = t - Cn a and N a = £ a + Cn a , (2.3) 

are respectively tangent and normal to the horizon. Note that V ■ V = — N ■ N = 2C. Hayward 
[4] has shown that if the null energy condition holds, then C > on a FOTH. Thus, the horizon 
must be either space-like or null, being null iff the shear as well as T ab £ a £ b both vanish 
across H. Intuitively, H is space-like in the dynamical regime where gravitational radiation and 
matter are pouring into it and is null when it reaches equilibrium. 

The second law of trapping horizon mechanics follows quite easily. If yjq is the area element 
corresponding to the metric q a b on the cross-sections, then 

Cv^q = -C9 {n) ^. (2.4) 

By definition 9i n ) is negative and we have just seen that C is non-negative, so we obtain the 
local second law: If the null energy condition holds, then the area element y/q of a FOTH is 
non-decreasing along the horizon. Integrating over St, the same law applies to the total area of 
the horizon sections. It is non-decreasing and remains constant if and only if the horizon is null. 
Dynamical Horizon: a smooth three-dimensional, space-like sub-manifold H of space-time 
is a dynamical horizon (DH) if it can be foliated by closed space-like two-manifolds St, with 
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future-directed null normals I and n such that (i) on each leaf the expansion 9^ of one null 
normal l a vanishes, (ii) the expansion 9r n \ of the other null normal n is negative. T 
Like FOTHs, a DH is a space-time notion defined quasi-locally, it is not relative to a space-like 
hypersurface, it does not refer to oo, it is not teleological. A space- like FOTH is a DH on 
which C n 9(c\ < 0; a DH which is also a FOTH will be called a space-like future outer horizon 
(SFOTH). The DH cannot describe equilibrium black holes since it is space- like by definition, 
but is better suited to describe how a BH grow in general relativity. Suitable analogues of the 
laws of black hole mechanics hold for both FOTHs and DHs. Our main interest in the following 
will be precisely for these local horizons, but for the time being we continue our description. 
Perfect and Non-Expanding Horizons: a perfect horizon is a smooth three-dimensional null 
sub-manifold H of space-time with null normal l a such that 6^ = on H and which intersect 
space-like hypersurfaces in compact sets. T 

If in the last clause H is topologically 1R x S 2 and moreover the stress tensor T a b is such that 
—T££ b is future causal for any future directed null normal £ a , then H is called a non-expanding 
horizon. T 

If X, Y are tangent to a non-expanding horizon we can decompose the covariant derivative 

V X Y = D X Y + N(X, Y)l + L(X, Y)n 

where Dx is the projection of the vector VjF onto the spheres St in H. If X is tangent to the 
spheres then Dx is the covariant derivative of the induced metric q a b, and if X is tangent to H 
one may regard the operator Vx = Dx + N(X, ■)£, acting on vector fields, as a connection on 
H. If this connection is "time independent" then the geometry of H is time independent too 
and we have Ashtekar et al. notion of an horizon in isolation. 

Isolated Horizon: a non-expanding horizon with null normal l a such that Vx] = along 
H.W 

These horizons are intended to model BHs that are themselves in equilibrium but possibly in a 
dynamical space-time. For a detailed description of their mathematical properties we refer the 
reader to Ashtekar-Krishnan's review [8]. 

The new horizons just introduced all have their own dynamics governed by Einstein eq.s. There 
exist for them existence and uniqueness theorems [17], formulation of first and second laws 
[8, 4, 18] and even a "membrane paradigm" analogy. In particular, they carry a momentum 
density which obey a Navier-Stokes like equation generalizing the classical Damour's equations 
of EHs, except that the bulk viscosity Qfoth = 1/167T > [19, 20]. We think Iver would be 
amused by that. 

The new horizons are also all space-like or null, hence it remains to see what is the role they 
play in the problem of black hole quantum evaporation. In this connection the following notion 
can be useful. 

Time-like Dynamical Horizon: a smooth three-dimensional, time-like sub-manifold H of 
space-time is a time-like dynamical horizon (TDH) if it can be foliated by closed space-like two- 
manifolds St, with future-directed null normals i and n such that (i) on each leaf the expansion 
6^ of one null normal i a vanishes, (ii) the expansion 6^ of the other null normal n is strictly 
negative. T 

Surface Gravities 

The surface gravity associated to an event horizon is a well known concept in black hole physics 
whose importance can be hardly overestimated. Surprisingly, a number of inequivalent defini- 
tions beyond the historical one appeared recently (over the last 15 years or so) in the field with 
various underlying motivations. We have collected the following (we rely on the nice review of 
Nielsen and Yoon [28]): 
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1. The historical Killing surface gravity (Bardeen et al. [21], textbooks) 

2. Hayward's first definition [4] 

3. Mukohyama-Hayward's definition [22] 

4. Booth-Fairhurst surface gravity for the evolving horizons [16] 

5. The effective surface gravity appearing in Ashtekar-Krishnan [8] 

6. The Fodor et al. definition for dynamical spherically symmetric space-times [23] 

7. The Visser [24] and Nielsen- Visser [25] surface gravity 

8. Hayward's definition [26] using Kodama's theory [27]. 

We will not spend much time on the various definitions and their motivations except for the last 
item, which is what the tunneling approach leads to, among other things. 

1. The Killing surface gravity is related to the fact that the integral curves of a Killing vector 
are not affinely parametrized geodesies on the Killing horizon H. Hence 

K a V b K a * KK a 

defines the Killing surface gravity n on H, where = means evaluation on the horizon. The Killing 
field is supposed to be normalized at infinity by K 2 = — 1. The definition can be extended to 
EHs that are not Killing horizons, by replacing K with the null generator I of the horizon. 
However there is no preferred normalization in this case, and this is one reason of the debating 
question regarding the value of the surface gravity in dynamical situations. 

2. Hayward's first definition was motivated by the desire to get a proof of the first law for THs. 
It is defined without appeal to inaffinity of null geodesies as 

^ = \^j-n a V a e {£) (2.5) 

and is independent on the parametrization of i a integral curves, since the evaluation is on a 
marginal outer surface where n ■ t = —1 and 9i(\ = 0. 

3. We leave apart the Mukohyama-Hayward and the Booth-Fairhurst definitions (see 4.) as 
they are somewhat more technical and complicated than it is necessary, so we refer the reader 
to the original papers. 

5. Given a weakly isolated horizon H, Ashtekar and Krishnan showed that for any vector field 
t a along H with respect to which energy fluxes across H are defined, there is an area balance 
law that takes the form 

5E l = —^—5 As + work terms 
8itG 

with an effective surface gravity given by 

1 dr 
K ~ 2RdR 

R is the areal radius of the marginally trapped surfaces, i.e. As = 4-7ri? 2 , the function r is related 
to a choice of a lapse function and finally E l is the energy associated with the evolution vector 
field t a . For a spherically symmetric DH a natural choice would be r = R so R = 1/2R, just the 
result for a Schwarschild BH. To illustrate the naturalness of this definition, consider a slowly 
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changing spherically symmetric BH with mass M(v), where v is a time coordinate. Defining the 
horizon radius at each time by R = 2M(v) and As = 4itR 2 , we can differentiate M to obtain 



which is the usual area balance law with surface gravity R = 1/2R = 1/4M. Consider, however, 
the more general possibility where the horizon is at R = 2M(v,R), as it happens for example 
in the Vaidya-Bardeen metric. The same computation leads to 

M = \ 2M'\ — =► « - 7TF f 1 - 2M ') ( 2 - 6 ) 



2R V ) 8vr AM 

a prime denoting the radial derivative. Thus naturalness is not a decisive criterion in this case. 

6. The Fodor et al. definition looks like the Killing form of surface gravity in that nl b = £ a V a £ b , 
where now £ a is an outgoing null vector orthogonal to a trapped or marginally trapped surface. 
This is because, as a rule, such null vectors are not affinely parametrized, although they can 
always be so parametrized that k = 0. So one needs to fix the parametrization: Fodor et al. 
choose 

K = -n a £ b V b £ a 

with n a affinely parametrized and normalized to n ■ t = — 1 at space- like infinity, t a being 
the asymptotic Killing field. Note that this definition is non local but looks like as a natural 
generalization of the Killing surface gravity. 

7. We postpone the discussion of the Visser and Visser-Nielsen surface gravity to the next 
section. 

8. Finally we have a local geometrical definition of this quantity for the trapping horizon of a 
spherically symmetric black hole as [26] follows. One can introduce local null coordinates x ± 
in a tubular neighborhood of a FOTH such that n = —g^ d- and £ = d + : then (shortening 

0{t) = #+> #(n) = 

1 

Later we will show that this k fixes the expansion of the metric near the trapping outer horizon 
along a future null direction. The definition may look somewhat artificial, but in fact it is 
very natural and connected directly with what is known for the stationary black holes. To see 
this one notes, following Kodama [27], that any spherically symmetric metric admits a unique 
(up to normalization) vector field K a such that K a G a b is divergence free, where G a b is the 
Einstein tensor; for instance, using the double-null form, one finds K = —g^ (d+rd- — <9_r<9+). 
The defining property of K shows that it is a natural generalization of the time translation 
Killing field of a static black hole. Moreover, by Einstein equations K a T ab will be conserved 
so for such metrics there exists a natural localizable energy flux and its conservation law. Now 
consider the expression K a V\pK a y. it is not hard to see that on H it is proportional to K^. 
The proportionality factor, a function in fact, is the surface gravity: K a V^K a ] = —kK^. For a 
Killing vector field VbK a is anti-symmetric so the definition reduces to the usual one. 

3 Two examples: Vaidya and McVittie's metrics 

We consider first spherically symmetric spacetimes which outside the horizon (if there is one) 
are described by a metric of the form 

ds 1 = -e 2 ^ r ^A(r, v)dv 2 + 2e*^Wr + r 2 dS 2 . (3.1) 



K=~{g+-d-0+),. =a (2.7) 
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where the coordinate r is the areal radius commonly used in relation to spherical symmetry 
and v is intended to be an advanced null coordinate. In an asymptotically flat context one can 
always write (we use geometrized units in which the Newton constant G = 1) 



A(r, v) = 1 — 2m(r, v)/r 



(3.2) 



which defines the active mass. This metric was first proposed by Vaidya [9], and studied in 
an interesting paper during the classical era of black hole physics by Lindquist et al [29]. It 
has been generalized to Einstein-Maxwell systems and de Sitter space by Bonnor- Vaidya and 
Mallet, respectively [30]. It was then extensively used by Bardeen [10] and York [11] in their 
semi-classical analysis of the back reaction problem. We will call it the Vaidya-Bardeen metric. 
A cosmological constant can be introduced by setting 



A(r, v) = 1 - 2m(r, v) jr - r 2 /L 2 



(3.3) 



where L~ 2 oc A. If one wishes the metric can also be written in double-null form. In the (v, r)- 
plane one can introduce null coordinates such that the dynamical Vaidya-Bardeen space-time 
may be written as 



ds 2 = -2f(x + ,x~)dx + dx- + r 2 (x + ,x~)dS 2 D _ 2 , 



(3.4) 



for some differentiable function /. The remaining angular coordinates contained in dS 2 do not 
play any essential role. In the following we shall use both forms of the metric, depending on 
computational convenience. The field equations of the Vaidya-Bardeen metric are of interest. 
They read 



dv dr 



-A-Kr 2 T v v , ^ = 47rre*T^ 
or 



(3.5) 



The stress tensor can be written as 



Tab = 



m m 



Airr 2 



2irr 2 



(3.6) 



If m only depends on v it describes a null fluid and obeys the dominant energy condition if 
m > 0. 

The second example we are interested in is the McVittie solution [12] for a point mass in a 
Friedmann-Robertson- Walker flat cosmology. In D-dimensional spacetime in isotropic spatial 
coordinates it is given by [31] 



ds 2 



-A(p, t)dt 2 + B(p, t) (dp 2 + p 2 dS 2 D ^ 2 ) 



with 



A(p,t) 



a(t)p 



D-3 



1 + 



(o(?)p) 



D-3 



B(p,t) =a{tf 



1 - 



m 
a(t)p 



D-3 



2/(D-3) 



(3.7) 



(3.8) 



When the mass parameter m = 0, it reduces to a spatially flat FRW solution with scale factor 
a(t); when a(t) = 1 it reduces to the Schwarzschild metric with mass m. In four dimensions 
this solution has had a strong impact on the general problem of matching the Schwarzschild 
solution with cosmology, a problem faced also by Einstein and Dirac. Besides McVittie, it has 
been extensively studied by Nolan in a series of papers [32]. To put the metric in the general 
form of Kodama theory, we use what may be called the Nolan gauge, in which the metric reads 



ds 2 = -(A s - H 2 {t)r 2 ) dt 2 + A~ l dr 2 - 2A-^ 2 H(t)r drdt + r 2 dS 2 D _ 2 



(3.9) 
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where H(t) = a/a is the Hubble parameter and, for example, in the charged 4-dimensional case, 

A s = 1 - 2m/r + q 2 /r 2 (3.10) 

or in D dimension A s = 1 — 2m/r D ~ 3 + q 2 /r 2D ~®. In passing to the Nolan gauge a choice 
of sign in the cross term drdt has been done, corresponding to an expanding universe; the 
transformation H(t) — > —H(t) changes this into a contracting one. In the following we shall 
consider D = 4 and q = 0; then the Einstein- Friedmann equations read 

3H 2 = 8tt P , 2Aj 1/2 H(t) + 3H 2 = -8tt P . (3.11) 

It follows that A s = 0, or r = 2m, is a curvature singularity. In fact, it plays the role that r = 
has in FRW models, namely it is a big bang singularity. When H = one has the Schwarzschild 
solution. Note how the term H 2 r 2 in the metric strongly resembles a varying cosmological 
constant; in fact for H a constant, it reduces to the Schwarzschild-de Sitter solution in Painleve 
coordinates. As we will see, the McVittie solution possesses in general both apparent and 
trapping horizons, and the spacetime is dynamical. However, it is really not a dynamical black 
hole in the sense we used it above, since the mass parameter is strictly constant: for this reason 
we called it a fake dynamical BH. This observation prompts one immediately for an obvious 
extension of the solution: to replace the mass parameter by a function of time and radius, but 
this will not be pursued here. 

The study of black holes requires also a notion of energy; the natural choice would be to use 
the charge associated to Kodama conservation law, but this turns out to be the Misner-Sharp 
energy, which for a sphere with areal radius r is the same as the Hawking mass [33], given by 
E = r(l — 2~ 1 r 2 g H 9 + 9-)/2. Using the metric (3.1) an equivalent expression is 

sTdnrdvT = 1 - 2E/r (3.12) 

In this form it is clearly a generalization of the Schwarzschild mass. As we said, E is just the 
charge associated to Kodama conservation law; as showed by Hayward [34] , in vacuum E is also 
the Schwarzschild energy, at null infinity it is the Bondi-Sachs energy and at spatial infinity it 
reduces to the ADM mass. 

Let us apply this general theory to the two classes of dynamical BH we have considered. Using 
Eq. (3.1), we have 9^ = A(r,v)/2r . The condition 9^ = leads to A(rh,v) = 0, which defines 
a curve = r^iy) giving the location of the horizon; it is easy to show that 9- < 0, hence 
the horizon is of the future type. Writing the solution in the Vaidya-Bardeen form, that is with 
A(r, v) = 1 — 2m(r, v) /r, the Misner-Sharp energy of the black hole is E = m(r^(v), v), and the 
horizon will be outer trapping if m (rh,v) < 1/2, a prime denoting the radial derivative. The 
geometrical surface gravity associated with the Vaidya-Bardeen dynamical horizon is 

k{v) - ^1 = m(r h v) _ m'(r h ,v) = 1 , _ a 
2 rt r h 4m V / 



' h 



the same as Eq. (2.6), where m = m(r/ l , v). We see the meaning of the outer trapping condition: 
it ensures the positivity of the surface gravity. 

As a comparison, Hayward's first definition would give k = vT— 2m' /4m, while Fodor et al. 
expression is 

2* 

k = — (l-2m') + * (3.14) 
4m 

The effective surface gravity of Ashtekar-Krishnan simply is k = l/4m, everything being evalu- 
ated on the horizon. Note that some of them are not correct for the Reissner-Nordstrom black 
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hole. We also note that k (3.13) is inequivalent to the Nielsen- Visser surface gravity, which in 
these coordinates takes the form 

k = — (1 - 2m - e~*m) (3.15) 

though they coincide in the static case. Also, both are inequivalent to the Visser surface gravity 
e^K, which was derived as a temperature by essentially the same tunneling method as discussed 
below, but in Painleve-Gullstrand coordinates. Part of the difference can be traced to a different 
choice of time. 

In the case of McVittie BHs, we obtain 

9± = ±(VA sT Hr)/2rf± 

where the functions f± are integrating factors determining null coordinates x^ such that dx^ = 
dt ± A s 1 ^ 2 dr \. One may compute from this the dual derivative fields d±. 

The future dynamical horizon defined by 8+ = 0, has a radius which is a root of the equation 
y/A s = Hrh , which in turn implies A s = H 2 r^. Hence the horizon radius is a function of time. 
The Misner-Sharp mass and the related surface gravity are 

E = m+ l -H(t) 2 rl (3.16) 

u\ 171 ij2 H E 3 2 H 
K{t) = ^ h - Hrh -2H = ^ h -2 Hrh -2H (3 ' 17) 

Note that E = In the static cases everything agrees with the standard results. The 

surface gravity has an interesting expression in terms of the sources of Einstein equations and 
the Misner-Sharp mass. Let T be the reduced trace of the stress tensor in the space normal to 
the sphere of symmetry, evaluated on the horizon H. For the Vaidya-Bardeen metric it is, by 
Einstein's equations (3.5), 

f = T v +T r = — — 

2irr h dr \ r =r h 

For the McVittie's solution, this time by Friedmann's equations (3.11) one has 
f = -p + p = -±(3H 2 ■ H ^ 



An \ Hr h 



We have then the mass formula 

kAh =E + 2nrlf (3.18) 



47T 

where Ah = Anr\. It is worth mentioning the pure FRW case, i.e. A s = 1, for which n(t) = 
- (H(t) + H/2H) . One can easily see that (3.18) is fully equivalent to Friedmann's equation. 
We feel that these expressions for the surface gravity are non trivial and display deep connections 
with the emission process. Indeed it is the non vanishing of k that is connected with the 
imaginary part of the action of a massless particle, as we are going to show in the next section. 
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4 Tunneling and instability 



The essential property of the tunneling method is that the action I of an outgoing massless 
particle emitted from the horizon has an imaginary part which for stationary black holes is 
3/ = 7tk~ 1 E, where E is the Killing energy and k the horizon surface gravity. The imaginary 
part is obtained by means of Feynman ie-prescription, as explained in [35, 36]. As a result 
the particle production rate reads V = exp(— 29/) = exp(— 2itk~ 1 E) . One then recognizes the 
Boltzmann factor, from which one deduces the well-known temperature T# = k/2tt. Moreover, 
an explicit expression for k is actually obtained in terms of radial derivatives of the metric on 
the horizon. 

Let us consider now the case of a dynamical black hole in the double-null form [37]. We have 
for a massless particle along a radial geodesic the Hamilton- Jacobi equation d + Id-I = . Since 
the particle is outgoing d-I is not vanishing, and we arrive at the simpler condition d+I = 0. 
First, let us apply this condition to the Vaidya-Bardeen BH. One has then 

2e-* ( - r ^d v I + A(r,v)d r I = 0. (4.1) 

Since the particle will move along a future null geodesic, to pick the imaginary part we expand 
the metric along a future null direction starting from an arbitrary event (r^(uo), vo) on the 
horizon, i.e. A(rh{vo), vo) = 0. Thus, shortening r^(vo) = ro, we have A(r,v) = d r A(ro,vo)Ar + 
d v A(ro, vq)Av + • • • = 2K(vo)(r — ro) + . . . , since along a null direction at the horizon Av = 0, 
according to the metric (3.1); here k(vq) is the surface gravity, Eq. (3.13). From (4.1) and the 
expansion, d r I has a simple pole at the event (ro,uo); as a consequence 

3/ = 9 a [* t „ 2e T' 8tJ T-^- («) 

J J A (r , v )(r - r - 10) k(v ) 

where u(vo) = e~^^ r °' v °^d v I , is to be identified with the energy of the particle at the time 
vq. Note that the Vaidya-Bardeen metric has a sort of gauge invariance due to conformal 
reparametrizations of the null coordinate v: the map v — > v(v), *$>(v,r) — > *S>(v,r) + ln(dv/dv) 
leaves the metric invariant, and the energy is gauge invariant too. Thus we see that the Hayward- 
Kodama surface gravity appears to be relevant to the process of particles emission. The emission 
probability, T = exp(— 2ttoj(v)/k(v)), has the form of a Boltzmann factor, suggesting a locally 
thermal spectrum. 

For the McVittie's BH, the situation is similar. In fact, the condition d+I = becomes 
d r l = -F{r,ty l d t I 

where 

F(r, t) = y/Ajr) [^AJrj - rH(tj) 

As before, we pick the imaginary part by expanding this function at the horizon along a future 
null direction, using the fact that for two neighbouring events on a null direction in the metric 
(3.9), one has t - t = (2H^y 1 (r - r ), where H = H(t ). We find the result 

F(r, t) = (± A' s (r ) - r F 2 - ^ (r - r ) + • • • = «(<d)(r - r Q ) . . . (4.3) 

where this time ro = r^to). From this equation we see that d r I has a simple pole at the 
horizon; hence, making use again of Feynman ie-prescription, one finds 3/ = 7TK(to) _1 w(to), 
where u(t) = dfl is again the energy at time t, in complete agreement with the geometric 
evaluation of the previous section. Obviously, if k vanishes on the horizon there is no simple 
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pole and the black hole should be stable 2 . The kind of instability producing the Hawking flux 
for stationary black holes evidently persists in the dynamical arena, and so long as the evolution 
is sufficiently slow the black hole seems to evaporate thermally. Note that the imaginary part, 
that is the instability, is attached to the horizon all the time, confirming the Fredenhagen-Haag 
suggestion quoted in the introduction. It is worth mentioning the role of k in the analogue of 
the first law for dynamical black holes (contributions to this problem for Vaidya black holes 
were given in [39]). Using the formulas of the projected stress tensor T given above, and the 
expression of the Misner-Sharp energy, one obtains the differential law 

provided all quantities were computed on the trapping horizon. Here Ah = ^kt h is the horizon 
area and Vh = 47rr^/3 is a formal horizon volume. If one interprets the "d" operator as a 
derivative along the future null direction one gets Hayward's form of the first law. But one can 
also interpret the differential operation more abstractly, as referring to an ensemble. Indeed, 
to obtain Eq. (4.4) it is not necessary to specify the meaning of the "d". It is to be noted 
that the same law can be proved with other, inequivalent definitions of the surface gravity, even 
maintaining the same meaning for the energy. Thus other considerations are needed to identify 
one: the tunneling method has made one choice. 

As thoroughly discussed in Hayward et al. [40], Eq. (3.1) is actually the most general form of a 
spherically symmetric metric, so the above calculations works throughout. Of course k > if 
the trapping horizon is of the outer type. Thus the method has derived a positive temperature 
if and only if there is a future outer trapping horizon. 



Extremal limit 



We discuss only an example, the charged stringy black hole, which represents a non-vacuum 
solution of Einstein-Maxwell dilaton gravity in the string frame [41, 42]: 

&2 - w+ (1 - a/r Ta-w -(^)'" 2 (45> 

where a > b > 0. The horizon radius is r = a. 

For this example, the extremal limit as defined by global structure is b — > a. The Killing surface 
gravity Kqo = l/2a does not vanish in this limit. Garfinkle et al. [42] noted this as puzzling, 
since extremal black holes are expected to be zero-temperature objects. 
Remarkably, the geometrical surface gravity (3.13) 

(4,) 

vanishes in the extremal limit. Thus the gravitational dressing effect lowers the temperature to 
its theoretically expected value. 

We conjecture that this is true in general. Indeed, past experience with extremal black holes 
showed that the horizon of these objects is not only a zero but also a minimum of the expansion 
9 + = d+A/A of the radially outgoing null geodesies, #_|_ becoming positive again on crossing the 
horizon. Thus d-9+ = should be the appropriate definition of an extremal black hole. Since 
k = —e~ 2,p d-d + r, this is equivalent to k = 0. 



2 However, charged extremal black holes can radiate [38]. 
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